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ABSTRACT 

On the basis of a large scale 'adiabatic', namely non-radiative and non-dissipative, 
cosmological smooth particle hydrodynamic simulation we compare the entropy pro- 
files of the gas and the dark matter (DM) in galaxy clusters. We employ the quantity 

K g = 3ksT g pg 2 ^ 3 /(/im p ) = <jgp g 2 ^ 3 as measure for the entropy of the intra-cluster 

gas. By analogy the DM entropy is defined as Xdm = o'dmPdm' 3 (°dm i s tne 3D veloc- 
ity dispersion of the DM) . The DM entropy is related to the DM phase space density 

by Kdm °c Qdm ■ In accord with other studies the radial DM phase space density pro- 
file follows a power law behaviour, Qdm oc r~ 182 , which corresponds to 1£dm oc r 121 . 
The simulated intra-cluster gas has a flat entropy core within (0.8 ± 0.4)i? s , where R s 
is the NFW scale radius. The outer profile follows the DM behaviour, K g oc r 121 , in 
close agreement with X-ray observations. Upon scaling the DM and gas densities by 
their mean cosmological values we find that outside the entropy core a constant ratio 
of K g /Kr)M = 0.71 ±0.18 prevails. By extending the definition of the gas temperature 
to include also the bulk kinetic energy the ratio of the DM and gas extended entropy 
is found to be unity for r > 0.8i? s . The constant ratio of the gas thermal entropy 
to that of the DM implies that observations of the intra-cluster gas can provide an 
almost direct probe of the DM. 
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1 INTRODUCTION 

Clusters of galaxies represent a composite of baryonic and 
dark matter (DM). The gas is commonly described by 
its thermodynamic variables, such as density, temperature, 
pressure and also entropy. The entropy provides a unique 
gauge for monitoring the thermal history of the gas. In the 
standard model of structure formation the primordial gas is 
isentropic. Ignoring radiative and dissipative processes and 
other non-gravitational sources or sinks of energy the only 
mechanism by which the gas entropy can change is via shock 
waves. According to the paradigm of hierarchical structure 
formation clusters form by merging and accretion of smaller 
building blocks, which usually contain baryons and DM. The 
collapse and merging of substructures lead to emergence of 
shock waves in the forming intra-cluster gas (ICM) and it 
ceases to be isentropic. These shocks can be seen close to 
the virial radius when the ga s is accreted onto the cluste r 
(see e.g. iBorgani et al.ll2005l: I Void 120051; IVoit et al.ll2005l ). 
However, the study of Pfrommer et alj (|2006l ) shows that a 



non negligible amount of energy is dissipated by shocks fur- 
ther inside the cluster potential well. In any case the result- 
ing entropy structure reflects the contribution of irreversible 
processes in the formation of clusters. The main problem ad- 
dressed here is to what extent the gas entropy resembles the 
structure of the DM halo and whether the classical thermo- 
dynamic description of the gas can be extended and applied 
to the DM. 

The entropy of the ICM can be deduced from X-ray ob- 
servations. As has become customary in this field, we refer 
to if oc T g pg 2 ^ s as the 'entropy', where T g is temperature 
and p g is the density of the gas. More precisely, we define 
K g = 3ksTgPg 2 ' 3 /(/im p ) = o-gPg 2 ^' i , where fes is the Boltz- 
mann constant and <r g is the internal 3D velocity dispersion 
of the gas. It is the latter expression upon which the direct 
comparison between the gas and the dark matter is based. 
Using XMM-New ton observations of 10 relaxed clusters 
iPratt et al.l (|2006h find the radi al profile of K follows a power 
law of r 13 with f3 « 1.1 (see also iDonahue et a"fll2006l h Apart 
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from the central regions the slope obtained i n adiabati c sim- 
ulations agree with the observed values (e.g. Voit 2005). The 
observations, however, show an excess in entropy which is 
thought to originate from non gravitational processes (e.g. 
iFinoguenov et al.ll20oj ; iBorgani et al.ll2005h . Note, the ac- 
tual entropy of an ideal gas is S = \n(K 3 ^ 2 ) + const., but 
here we follow the common convention and associate the 
entropy with K. 

The concept of entropy has also been applied to col- 
lisionless self-gra vitating systems in its Boltzmann-Gibbs 
formu lation (e.g. iLvnden-Belll 1 19671 ; iLvnden-Bell fc Wood! 
1968) and in a n analogy with the thermodynamic entropy of 
a perfect gas (IWhite fc Naravanlll987l ; iNavarro et~ai1ll995l ; 
lEke et all Il99a ). It is the latter approach that we follow 
here. In analogy with classical thermodynamics and the 
ideal gas equation of state we write the local specific en- 
tropy 5dm = 1ii(od M /pdm), where ctdm is the DM velocity 
dispersion and an additiv e constant i s negl ected. A similar- 
approach is presented by lEke et al.l (|l998l ). Following the 
convention used in the studies of clusters of galaxies we de- 
fine -Kdm = °"dmPdm 3 - Note that the DM specific entropy 
is the log of the inverse of the so-called phase space den- 
sity Qd m = PDingnM- Now , the a nalysis of simulated DM 
halos by ITavlor fc Navarrd (|200ll ') has shown that the ra- 
dial phase space density profile closely follows a power law, 



Qdm(t) <x rj 



Using a set of high-resolution N-body 



simulations lAscasibar et al.l ([2004) confirm this power-law 
behaviour for an extend range of haloe masses, in addition 
they find coincidence of the zero-points if the profiles are 
conveniently rescaled by virial quantitie s . Cha nging vari- 
ables from Q to K the lTavlor fc Navarrd (|200lh power law 
is given by 



^DM(r) = QdmV) « r 



(1) 



The close simila r ity o f the (appropriately converted) 
ITavlor fc Navarrol (|200ll) slope an d the slope derived from 
X-ray observations by Pratt et al.l (|2006T ) has motivated us 
to compare the DM and gaseous entropy profiles of clus- 
ters of galaxies extracted from a large scale 'adiabatic' 
Smoothed Particle Hydrodynamics (SPH) simulation. In 
this context adiabatic is defined as non-radiative with no 
non-gravitational sources or sinks of energy. In such kind 
of simulations the entropy of each particle remains constant 
as lo ng as it does not pass through a shock front (|Springell 
2005). The only dissipative process is introduced by the nu- 
merical viscosity which damps unphysical oscillations be- 
hind shock fronts and thus transfers kinetic into thermal 
energy in accordance with the Rankine-Hugoniot jump con- 
ditions. 

This paper is organised as follows. Sj2] introduces the 
simulation and the cluster ensemble and Sj3] compares the 
entropy profiles of the gaseous and the DM component of 
the clusters. i|4] focuses on the emerging gas entropy cores 
and discusses their physical properties. And a discussion is 
given in SJS] 



2 SIMULATION 

This simulation, dubbed Th e MareNostrum Uni verse 
and described in detail in Gottlob er et al.l ((2006), was 
performed with the entropy conserving Gadget2 code 



l|Springeli r2005V It followed the non linear evolution of the 
cosmic density fields of the gas and the DM from z — 40 
to the present epoch (z — 0) within a comoving cube of 
500 h _1 Mpc edges. The spatially flat concordance model em- 
ploying the following parameters was adopted: O m = 0.3, 
Qa = 0.7, Q b = 0.045, the Hubble parameter h = 0.7, 
the slope of the power spectrum n = 1 and the normali- 
sation erg = 0.9. Both components, the gas and the DM, 
are resolved by 1024 3 particles, which results in a mass of 
8.3 x 10 9 h _1 M for the DM particles and 1.5 x 10 9 h _1 M 
for the gas particles. The gas dynamical equations were 
solved by SPH without including dissipative or radiative 
processes or star formation, i.e. the gas was treated adia- 
batically in the above discussed sense. The SPH smoothing 
length has been set to the 40 th nearest neighbour of each 
SPH particle and a gravitational force resolution was cho- 
sen equivalent to a Plummer softening of 15h _1 kpc (comov- 
ing). The SPH softening was not allowed to fall below the 
gravitational softening length. 

The cluster ensemble comprises 96 massive clusters se- 
lected from the simulation volume. The clusters were found 
by a friends of friends algorithm (FoF) with a linking length 
of 0.17 times the mean particle separation applied to the 
DM component only. This corresponds roughly to the virial 
overdensity. In a second step substructures within the clus- 
ter volumes were identified by reducing the linking length of 
the FoF group finder by a factor of 4 which defines substruc- 
tures at roughly 64 times the virial overdensity. The choice 
of this particular linking length is somewhat arbitrary, how- 
ever it has been found to adequately identify relevant sub- 
structure, such as double or triple core systems. The cluster 
centre was associated with the centre of the most massive 
substructure of the cluster. Subsequently the virial radius 
-Rvir determined as the cluster centric radius which encloses 
334 times the mean cosmic density and the virial mass were 
computed as the sum of the gas and the dark matter mass 
within that radius. The virial masses lie in between 3.8 x 10 14 
and 2.5 x 10 15 h _1 M Q corresponding to ~ 5 x 10 4 - 2 x 10 5 
particles per species. In total 370 clusters with masses above 
3.8 x lO 14 h _1 M are found within the simulation volume 
and for a even lower limit of 10 14 h _1 M we obtain 3708 ob- 
jects. The selection process for the clusters in combination 
with the steep mass function causes our sample to thin out 
at the low mass end, but for high masses (> 5 x 10 14 h _1 M ) 
it can assumed to be complete. To assure the best possible 
determination of the central properties of the cluster the 
centre was redefined as the associated peak of the density 
field smoothed by a 100h _1 kpc sphere. On average this re- 
adjustment changed the location of the centres by only a 
few h _1 kpc. 



3 ENTROPY PROFILES OF THE GAS AND 
THE DARK MATTER 

The aim of this section is to demonstrate the similarities in 
the entropy profiles of the dark and the gaseous matter. As 
motivated above we will refer to the entropy of the gas and 
the DM as 



Ki 



2/3 



(2) 
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Figure 1. Radial profiles for three arbitrarily selected clusters. The uppermost panels display the density profiles of the DM (solid line) 
and the gas (dashed line) scaled by their mean cosmic values. Here c indicates the DM concentration derived from a NFW-fit (dotted 
line). The second row shows the mean radial velocities for the dark matter (solid line) and the gas (dashed line). The third row displays 
the velocity dispersions a. The DM dispersion is given by the solid line. The turbulent dispersion of the gas is shown as dashed line. 
The thick dashed-dotted line depicts the thermal velocity dispersion of the gas and the thin dashed line displays the extended one. The 
lowest panels show the DM and gas entropy. The same line styles and units are used as in the panels above. 



were p is the density, a is the three dimensional velocity dis- 
persion and i stands for either DM (dark matter) or g (gas). 
The velocity dispersion of the dark matter within a radial 
bin is computed as the rms velocity of the particles within 
that bin. By their nature, the gas particles are characterised 
by two distinct properties related to their kinetic energy, 
namely, a proper velocity (same as for dark matter particles) 
and a temperature. The proper motions of the gas particles 
will be referred to as bulk motions and the corresponding 
velocity dispersion will be addressed as turbulent dispersion. 
Also the temperature of the gas will be expressed by means 
of the velocity dispersion that represents the thermal energy 
per unit mass, namely by a 2 = 3fcsT/(/im p ), where m p is 
the proton mass and [i — 0.588 for a completely ionised gas. 
The velocity dispersion of the gas derived from its tempera- 
ture will be referred to as thermal dispersion. The composite 
of the two gaseous dispersions is the extended dispersion. 

Before proceeding to show the average properties of the 
cluster ensemble we examine in some detail the structural 
and thermal properties of three typical clusters, selected 
with the intention to approximately cover the range of con- 
centrations apparent in the ensemble. Fig. [T] displays radial 
profiles, namely the density, radial bulk velocity, velocity 



dispersion and entropy profiles, of three arbitrarily chosen 
clusters of comparable masses (~ 10 15 h -1 Mo). The upper- 
most panels show the density profiles for the DM and the 
gas scaled by the corresponding cosmic mean density . The 
concentration parameter c = R v j T /R a specified in the up per 
right corners is based on a NFW-fit <|Navarro et al.ll 19971 ) to 
the DM component and the associated profile is displayed 
with the dotted lines. The gas velocity dispersion and en- 
tropy is decomposed into the thermal and kinetic compo- 
nents. In addition the extended dispersion and entropy of 
the gas is shown. 

The important features to notice in Fig. [T] are: (a) The 
excellent NFW fit to the DM density profile of the high 
concentration cluster, while the gas density exhibits a core 
structure; (b) The lack of significant radial streaming mo- 
tions of the high concentration cluster implies that it is a 
relaxed cluster while the others are still dynamically evolv- 
ing; (c) The DM velocity dispersion profile of all clusters 
show a 'temperature inversion', namely the velocity disper- 
sion increases with radius within « 0.1-Rvir- The gas thermal 
and extended velocity dispersions at the very centre (below 
« 0.05-Rvir) show the opposite behaviour, at these radii a 
negative temperature gradient of the gas is apparent in all 
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Figure 2. The upper two panel display the DM and the gas 
entropy profiles of the clusters within the ensemble. To reduce the 
scatter caused by the varying masses of the clusters the densities 
have been scaled by their mean values within the virial radius 
and the velocities have been scaled by the virial circular velocities 
(GMvir/iJvi,.) 1 / 2 . The lowest panel shows the mean DM and gas 
entropies by the thick solid and dashed-dotted lines. The thin 
dashed-dotted line is the mean profile of the extended entropy of 
the gas. The dotted line on top of the solid line represents the best 
fitting power law for the DM component with a logarithmic slope 
of /3 = 1.21. The corresponding exponent for the phase space 
density is a = —1.82. 



the clusters. Cluster c = 6.9 displays an almost isothermal 
core with the hint for a maximal temperature at the same 
radius where the DM dispersion peaks, (d) The DM entropy 
profile is close to a power law inde pendent of the dynam i- 
cal state of the cluster (compare to lAscasibar et all 12004 ). 
Within ~ 0.1i?vir the thermal and extended gas entropies 
show a core structure. Outside that region the extended gas 
entropy almost perfectly coincides with the DM entropy. 

Fig- El presents a direct comparison of the DM and 
the thermal and extended gas entropy profiles as defined in 
Eq. [2] for the 96 clusters in the sample. The DM entropy is 
displayed in the upper panel and the gas thermal entropy 
is given in the middle panel. The bottom panel exhibits the 
mean of the DM and the thermal and extended gas entropy 
profiles. The mean profiles reproduce the trends followed by 
the individual clusters of Fig. [T] The DM profile follows a 
power law of the form Kum tx r 13 , where (5 — 1.21. This 
logarithmic slope corresponds to a slope for the phase space 
density profile of a = —1.82. The mean gas thermal and 
extended profiles keep almost a constant displacement (in 
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Figure 3. Ratio of gas and DM entropies. The mean value of 
Kg/ K DM a t the radius indicated by short vertical line is dis- 
played by the dashed horizontal line. 

the log- log plane). The gas exhibits an entropy core of a 
radius R CO ic ~ 0.1i? v ir- Outside that core the gas entropy 
trails the DM power law. In particular, the extended entropy 
coincides with the DM entropy. The lower panel of Fig. [2] 
presents the main new result of the paper, namely that the 
gas extended entropy outside the core coincide with the DM 
entropy not only in slope but also in amplitude. 

The thermal entropy of the gas is more accessible to 
astronomical observations than the extended one. So, the 
ratio of the thermal to the DM entropy profile of all the 
clusters is presented Fig. [3] Outside the core that ratio is 
virtually constant and exhibits a small scatter, K s /Kdm = 
0.71 ±0.18. 



4 GAS CORE STRUCTURE 

Th e pre s ent re sults provide further support for the findings 
by IVoitl (|2005h that the emergence of the entropy core in 
adiabatic simulations is a real physical effect and is not a nu- 
merical artifact. Another indication for the reality of the en- 
tropy core comes from the compa rison of SPH with Eulerian 
grid codes. Ascasiba r et al.l (|2003h found good agreement be- 
tween the same cluster simulated by both techniques. In par- 
ticular an entrop y core emerged in both realisations. How- 
ever, as noted by IVoit et all |2005l) the c ore entropy levels 
can differ by almost a factor of two. Also, iLin et ah! (|2006f ) 
point out that the resolution of the SPH simulations plays 
an important role. Yet, the resolution for the ensemble pre- 
sented here is sufficient to exclude the pure artificial nature 
of the entropy cores. This statement is based in particular 
on the intermediate and low concentration systems. Here the 
radius of the core region is > 200 h~ kpc (assuming a virial 
radius of 1 h _1 Mpc and a concentration of 5) which can not 
be a caused by a lag of spatial resolution. It follows that the 
mere existence of entropy cores does not imply non-adiabatic 
physics such as feedback and preheating. 

To improve our understanding of that phenomenon we 
undertake a more detailed investigation of entropy core. As 
we are interested in the self-similar properties of the clusters 
all further plots of radial profiles will be scaled by the NFW 
scale radius R s , rather than by i? v ir- At the same time the 
density profiles are scaled by their mean densities within R s 
and the velocities are scaled by the circular velocity at R s , 
v c ,s = (GMs/Rs) 1 / 2 , where M s is the total mass within R s . 
That kind of scaling achieves a remarkable agreement among 



Entropy of gas and dark matter in galaxy clusters 5 




Figure 4. The upper two panel show the temperature profiles for 
the DM and the gas, respectively, of the clusters sample. The pro- 
files are plotted against the radius scaled by the NFW R a . The 
different line styles of the profiles are chosen according to the 
concentration parameter, where clusters are grouped into c < 5, 
5 < c < 6 and 6 < c. The bottom panel shows the mean pro- 
files, where the mean is taken over the groups with the above c 
range. Both the thermal and extended temperature profiles are 
presented. The DM 'temperature' (namely velocity dispersion) 
shows a clear temperature inversion within 0.8-R s . The gas pro- 
files is slightly monotonically decreasing within that inner region. 



the dark matter density profiles, implying that the scatter of 
the gas density profiles, in particular at the centre, is a real 
physical effect (Figs. U [S] and HJJ. To quantify the scatter 
in the density and velocity dispersion profiles we compute 
the rms value of (x — x)/x, where x indicates the values 
associated with the individual profiles and x is their mean 
value within a particular bin. As shown in Fig.|6]the scatter 
of the DM density lies between 0.2 and 0.1 and a minimum 
is located near R s . The central scatter p gas is 0.6 but it 
declines with radius and converges to the DM scatter at 
~ 0.8R S . The scatter in a dm is very similar to that of the 
DM density profile, ranging in between to ~ 0.2 to ~ 0.1. 
The scatter in the thermal velocities of the gas ranges from 
0.5 at the centre to 0.3. ft follows that there is a much larger 
variation among the different clusters in the gas density and 
temperature profiles than in the DM component. 

To explore the dependence of the core structure on the 
dynamical state of the cluster the ensemble has been subdi- 
vided into three concentration bins. The three bins, namely 
c<5, 5<c<6 and 6 < c comprise 39, 19 and 38 clus- 
ters, respectively. The intermediate bin should be considered 



Figure 5. The upper two panel show the scaled density pro- 
files for the DM and the gas, respectively, of the clusters sample. 
The profiles are plotted against the radius scaled by R a . The line 
styles follow the one of Fig. [4] The bottom panel displays the 
mean profiles, where the mean is taken for each concentration 
bin separately. 
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Figure 6. Radial dependence of the scatter for the density and 
velocity profiles of the gas and the DM. Displayed are the rms 
values of (x — x)/x, where x indicates the values associated with 
the individual profiles and x is their mean value within the bin. 
The thick solid and dashed lines show the scatter in the gas and 
DM densities (Fig. [5}. The thin solid and dashed lines depict the 
scatter in the gas and DM velocity dispersions (Fig. [4). 



as a kind of a buffer zone for the statistically more robust 
bins on the edges. The clusters with concentrations below 
5 are assumed to be on average dynamically young systems 
whereas the clusters with concentrations above 6 are ex- 
pected t o be the most relaxed systems within the ensemble 
see e.g. IWechsler et~alll2002i ). 
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Figure 7. The mean DM and gas thermal entropy profiles is 
presented for the different concentration parameter grouping. The 
DM entropy file follows the r oc r@ power law, independent of the 
value of c. 



Figure 8. Distribution of i? ma x/K s vs R B , where i? m ax is the 
radius at of the maximum of the DM 'temperature'. The distri- 
bution is characterised by i? ma x/J2a = 0.8 ± 0.4. 



Under the assumption that the ICM is in hydrostatic 
equilibrium the gas entropy excess within Reave implies that 
the gas is hotter and/or the gas density profile is flatter than 
the DM. This is confirmed by the simulation. The radial 
structure of the DM and gas temperature is given in Fig. [4] 
The figure nicely illustrates the temperature inversion of the 
DM temperature profile and, in contrast, the gentle increase 
of the gas temperature towards the centre. Fig. [5] presents 
the DM and gas density profiles and it clearly shows the 
flattening of the inner gas density profile. The gas entropy 
core is established by the flattening of the NFW density cusp 
and the absence of the temperature inversion. 

Next we focus on the impact of the concentration pa- 
rameter on the entropy profiles. Fig. [7J shows the mean 
DM and gas thermal entropy profiles and their dependence 
on the value of c. The DM is found to be independent of 
c and the gas shows a very marginal dependence. A much 
stronger variation of the core entropy level and the core sizes 
with concentrations is seen if i?„i,.-scaled quantities are used 
(not shown here). In that context it is worth noting, that 
the mean density profiles for the different concentrations do 
match very well apart from the inner 0.3R B where the low- 
est concentration profiles are slightly enhanced (Fig. [5}. In 
contrast the mean temperature profiles of the gas (Fig. [4} 
show a slight dependence on concentrations, namely the gas 
temperature rises with concentration. This trend is not re- 
flected in the behaviour of the DM velocity dispersion pro- 
files and therefore it may be real and not caused by the v C:3 
scaling. Therefore, higher concentrated clusters have slightly 
higher gas temperatures. Now, rising temperatures at equal 
densities result in higher entropies for higher concentrated 
clusters as seen in Fig. [7] However, it can be questioned 
whether the scaling of the velocities by the circular velocity 
at Rs (v c ,s) is the physically appropriate way of gauging the 
temperature profiles of different sized clusters. 

The core radius seems to coincide with the radius at 
which the DM temperature reaches its maximal value, de- 
fined as -Rmax- Thus it seems that the 'thermal', namely 
velocity dispersion of the DM, properties of the DM either 
set, or otherwise are related to, the extent of the entropy 
core. Fig. [8] presents a scatter plot of the ratio Rmax/Rs 
versus R B . There is no significant dependence of R ma x/Rs 
on R s or on c. Associating R mBX with the core radius we 
conclude that R covc /R s = 0.8 ± 0.4. 



5 DISCUSSION 

The main findings of the paper are: The DM phase space 
density power law of lTavlor fc Navarro! (|200lh is reproduced 
by the simulated clusters extracted from an adiabatic gas dy- 
namical simulation. The gas entropy p rofile is in accord with 
previous simulations l)Voit et al.ll2005h . In particular the en- 
tropy profile has an inner flat core and an outer profile of 
K g oc r 1 ' 21 is recovered. This stands in at least qualitative 
agreement with X-ray observations. The entropy core coin- 
cides with the region of the temperature inversion of the DM 
and its radius is characterised by R CO rc/Rs = 0.8±0.4. How- 
ever, entropy loss due to radiative processes, which is not 
included in our simulati ons, alters the thermo dynamics of 
the central ICM, see e.g. iKravtsov et"al] (120051). Therefore, 
it is n o contradiction to our findings that Piffaretti et al.l 
(2005) don't detect entropy cores within 13 nearby cool- 
ing flow clusters. The most important result of the paper 
is that outside the entropy core the gas and DM entropies 
follow one another very closely. Upon appropriate scaling of 
the gas and DM densities by their cosmological mean values 
we find that outside the entropy core a constant ration of 
K s / Kum = 0.71 ±0.18 prevails. By extending the definition 
of the gas entropy to include also the kinetic energy of the 
turbulent motion the above ratio becomes unity, and it fol- 
lows that in the outer part of the cluster the gas and DM 
entropies coincide and are virtually the same. 

The constant ratio of the gas to the DM entropies out- 
side the entropy core has interesting observational ramifica- 
tions. To the extent that non-adiabatic physics, in the sense 
described above, does not affect the gas structure outside the 
entropy core observations of the ICM can shed light on the 
DM in clusters. Observational determination of the entropy 
of the ICM can be easily translated into the entropy of the 
DM. Given the power law nature of the DM entropy, the ob- 
servationally determined outer DM entropy can be extrapo- 
lated inwards. This, together with independent mapping of 
the mass distribution can yield interesting insight into the 
distribution and dynamics of DM in clu sters of galaxies. In 
this context it is worth mentioning that llkebe et alj (|2004| ) 
described a formalism how the DM "temperature" can be 
derived from X-ray observations and applied this new tech- 
nique to XMM-Newton data of A1795. 

The analysis of § 0] shows that the entropy core radius, 
measured in R s units, is independent of the cluster NFW 
parameters. This seems to imply that the emergence of the 
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entropy core is related to processes occurring within 7? a . 
The common understanding of halo formation is that the 
structure within R s is determined exc lusively in the violent 
merger phases the cluster undergoes (jWechsler et al.l 120021 : 
iRomano-Diaz et al.ll2006h . The quiescent accretion phases 
that follows the violent mergers do not alter the inner struc- 
ture but leads to a continuous growth of virial radius. This 
suggests that the entropy core emerges as a result of pro- 
cesses occurring during the violent phases, in which the clus- 
ter is not in NFW equilibrium. Another alternative is that 
the gas component, unlike the DM, undergoes some secu- 
lar evolution of the inner region in the quiescent phases. It 
might also be the case that both processes play an impor- 
tant role. In any case it is clear that only shock waves can 
form the entropy core and the understanding of the shock 
pattern is critical here. 

N-body simulations show that clusters have experi- 
enced merging events in their past. The merging events 
of comparable sized substructures are violent enough to 
send bow shocks into the intra-substructure gas which 
subsequently propagate to the outskirts of the clusters 
dNagai fe Kravtsovll2003l ) . This shocks raise the entropy of 
the gas. IPfrommer et al.1 (|2006l ) develop a formalism for the 
identification and accurate estimation of the strength of 
structure formation shocks during cosmological SPH sim- 
ulations. They find that most of the energy is dissipated in 
weak internal shocks with Mach numbers M = 2 which are 
predominantly central flow shocks or merger shock waves 
traversing halo centres. The study by iNavarro et al.1 (|l995l ) 
utilises cosmological simulations to compare the evolution of 
DM and gas of clusters between z = 1 and z — 0. They find 
that the DM evolves at roughly constant phase space den- 
sity, whereas shocks increase the t emperature of the gas w ith 
relatively little change in density. lAscasibar et al.l (I2004T ) in- 
vestigating pure N-body simulations report the remarkable 
stability of the phase space den sity profile of the D M inde- 
pendent of the dynamical state. IPearce et al.1 (|l994h explore 
simple head-on mergers between model galaxy clusters con- 
taining a mixture of gas and dark matter, testing the cou- 
pling of the gas to the underlying collisionless material. Af- 
ter the first passage of the dark matter cores the gas forms 
a pancake-like structure about the centre of mass. While 
separating for a second time dark matter cores leave their 
gas behind. The central gas is no longer confined to a deep 
potential well so it expands. This pro cess has been nicel y 
confirmed by the recent observations of lClowe et alj (2006). 
They find an obvious displacement of the gas measured by 
X-ray observations compared to the potential revealed by 
weak lensing maps of the merging cluster system IE 0657- 
558. At the same time gravitational attraction due to the 
central gas exerts a drag force on the separating DM cluster 
cores, c ausing an energy transfe r from the DM to the gas 
(see also INavarro fc Whitell 19931). This is in agreement with 
recent investigations by iLin et ail (|2006l ) reporting a ~ 10 
per cent increase in the concentration of the cluster's DM 
component when adiabatic gas physics is included as com- 
pared to pure DM simulations. The core structure of the 
gas entropy, in contrast to the power-law profile of the DM, 
may therefore be established during merging events in con- 
sequence of the impenetrability of the gas as opposed to the 
collisionless nature of the DM. 

A different point of view might be provided by the fol- 



lowing speculations. The DM 'temperature' profile is char- 
acterised by a temperature inversion. Now, a 'normal' ther- 
modynamic system with a mechanism for heat transport 
and no energy sinks or sources cannot sustain such a struc- 
ture. Heat would flow inwards, leading to an increase of the 
entropy where the temperature inversion prevails. The in- 
crease in entropy amounts to an increase in temperature, 
and thereby decrease in density so as to keep the hydrostatic 
equilibrium. Now, let us conjecture that the collapse of the 
gas and DM results in a configuration where the gas follows 
the DM distribution. (This is supported by the self-similar 
solutions of Bertschinger 1985 and Chuzhoy & Nusser 2000). 
If the gas had a mechanism of transporting energy it would 
behave in the manner described above, namely the entropy 
will increase where the initial temperature inversion initially 
existed. Inspection of Fig. [T] proves that this might have ac- 
tually happened. The gas entropy departs from the DM and 
an entropy core replaces the temperature inversion. This is 
clearly manifested by the most relaxed cluster of Fig. [1] and 
is seen in other relaxed clusters (not shown here). So, what 
can provide the heat transport mechanism? Fig. [T] shows 
that outside the entropy core about (25 — 30) per cent of 
the non-gravitational energy is in the form of kinetic energy. 
Only a small fraction of this energy is in large scale bulk mo- 
tion and therefore most of it is in small scale motion, namely 
turbulence. The dissipation of these motions by weak shocks 
is expected to be efficient at the centre and this would lead 
to the increase in the gas entropy. In the language of one di- 
mensional spherically symmetric hydrodynamics one would 
attribute this phenomenon to turbulent viscosity and heat 
transport. This transport mechanism erases the temperature 
inversion and produces the entropy core. 

The application of the concept of maximum en- 
tropy as a mean of calculating the equilibrium configu- 
ration of s elf gravitating systems has proven somewhat 
futile (e.g. iLvnden-Belil Il967l: iLvnden-Bell fc Wood! Il968l ; 
IWhite fc Naravanl Il987l iLima Neto et al.lll999T ). In partic- 
ular, these attempts have not predicted the NFW density 
profile that faithfully reproduces the structure of DM ha- 
los. The problem that lies at the root of the shortcoming 
of the thermodynamic approach is the long range nature 
of the gravitatio nal force and the resulting n egative specific 
energy (see e.g. iLvnden- Bell fc Woodl H968). In spite of all 
that our results show that the concept of the DM entropy 
is closely related to the structure of DM halos. Moreover, 
the DM entropy is very closely related to the gas entropy. 
Beyond the entropy core of the gas the entropies of the two 
components are tracing one another. In fact, the extended 
gas entropy fully coincides with the DM entropy in that 
regime. It follows that the observed power law behaviour of 
the gas entropy provides an indirect support for the NFW 
structure of the DM. Now, this opens the door for a funda- 
mental question. The gas entropy results from the conver- 
sion of gravitational to kinetic energy and by the dissipation 
of the kinetic energy into internal one by means of shock 
waves. The power law regime seems to be dominated by the 
accretion shock that is predicted by self-sim ilar solutions 
l|Bertschingeijll985l ; IChuzhov fc Nussed I2OO0T ). In the DM 
sector the gravitational energy is converted to random mo- 
tions by means of phase mixing and the time variation of the 
gravitational potential. The latter process happens by vio- 
lent relaxation, whereby the gravitational potential changes 
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rapidly, or by means of a diabatic changes as envisaged in the 
secondary infall model (Gunn 1977 ; Fillmore fe Goldreichl 
1 19841 ; iBertschingerl 11983 : iHoffman fe Shahaml Il985l ). The 
physical processes that shape the entropy profile in the gas 
and the DM appear to be very different, yet the resulting 
structure seems to be identical. This phenomenological find- 
ing needs to be further investigated. 

Additional studies of the adiabatic formation of clus- 
ters of galaxies are needed to deepen our understanding of 
the role of entropy in the formation and evolution of clus- 
ters of galaxies. In particular, high resolution gas dynamical 
simulations need to be analysed with a special emphasise 
on the role of strong, large scale, shocks and of weak, small 
scale, ones. The relevance of these to the emergence of the 
entropy core needs to be further clarified. Finally, the sim- 
ulations should be expanded so as to include preheating, 
radiative cooling and feedback processes. 
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